Spin projected unrestricted Hartree— Fock ground states for harmonic quantum dots. 
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We report results for the ground state energies and wave functions obtained by projecting spatially 
unrestricted Ifartree Fock states to eigenstates of the total spin and the angular momentum for 
harmonic quantum dots with A*" < 12 interacting electrons including a magnetic field. The ground 
states with the correct spatial and spin symmetries have lower energies than those obtained by 
the unrestricted method. The chemical potential as a function of a perpendicular magnetic field is 
obtained. Signature of an intrinsic spin blockade effect is found. 

PACS numbers: 73.23.Hk, 73.63.Kv 
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I. INTRODUCTION 

Systems like atoms, metal clusters, trapped bosons and 
quantum dots show several universal features i For ex- 
ample, strongly interacting electrons in quantum dots 
arrange themselves in a rotating Wigner molecule.— Ro- 
tating boson molecules have been predicted to exist in ion 
traps.— Furthermore, symmetric potentials can induce a 
shell structure in atoms, metal clusterSf^ and quantum 
dots.i^ In the latter, signatures of shell structure have 
been experimentally probed,—"^ leading to Hund's rules 
for the total spin of the electron ground state. The spin 
in quantum dots'* also affects the electron transport. It 
can lead to spin blockade effects^ii^ and negative differ- 
ential conductance in nonlinear transport and 
it induces periodic modulations of the positions of the 
Coulomb peaks in the linear conductance as a function 
of an applied magnetic field. ^-'^■-'^^•^^■^Z' More recently, the 
effect of the spatial distribution of the spins on the Kondo 
phenomenon has been probedj ^^i^^'^°i^^ 

Electron and spin states of quantum dots have 
been theoretically studied with various techniques.— 
For small electron numbers N, exact diagonalization 
(ED),^>23^>25^>27^>29^^^^ configuration interac- 
tion (CI);^'^ and stochastic variational methods^ allow 
for determining ground and excited state energies and 
their quantum numbers with high accuracy. For larger 
N, the size of the many-body basis set increases exponen- 
tially. With presently available computational technol- 
ogy, reliably converged "exact" results can be obtained 
only for electron numbers up to TV « 8 electrons 1^ 

For < 13, iV = 16,24,48, Quantum Monte 
Carlo (QMC)^'^i29do.4i^d3^ methods have been 
used. They can provide accurate estimates for ground 
and excited states energies. With these techniques, the 
shell structure, Hund's rules, Wigner crystallization and 
the occurrence of "magic" angular momenta have been 
investigated^'^'^'^ii^'i^iil'i^ Most of the results for 
higher particle numbers have been restricted to zero mag- 
netic field. It is believed that QMC provides better es- 



timates for the energies of the ground states for larger 
electron numbers as compared to the "exact" methods. 

For larger N and/or in the presence of magnetic field 
B, methods like Hartree Fock (HF) 49^^^^^^ ^nd 
the density functional theory^^^i^i^i^iSiiS^ have been 
used. Generally, these seem to provide less accurate es- 
timates for the ground states which also can have un- 
physical broken symmetries due to incomplete ansatz 
wave functions. For instance, neglecting correlations, the 
straightforward HF method starts from a single Slater 
determinant as a variational many-body wave function 
which not necessarily is an eigenstate of the total spiujS 
Spatially unrestricted HF methods (UHF)— 1^ systemat- 
ically use symmetry breaking in order to obtain better 
estimates for the ground state energy. This may lead to 
wrong results for the total angular momentum and the to- 
tal spin. For instance, UHF calculations sometimes seem 
to fail predicting the total spin resulting from Hund's 
rule, in contradiction to the more accurate methods. Vi- 
olations of Hund's rules for relatively weak Coulomb in- 
teractions have been reported**^ for = 4, 8, 9. 

Projection techniques, pioneered in the 60th of the last 



century 



64.65.66 



can be applied for introducing the correct 



spatial symmetries. In quantum dots, they have been 
used for obtaining wave functions corresponding to spe- 
cific angular momentaj ^^i^^'^^i^"'^^ Recently, the random 
phase approximation has been used to restore the rota- 
tional symmetry of wave functions obtained by UHFi^ 
Restoring the spin symmetry has received much less 
attention, and has been used only for very few (up to 
iV = 3) electrons. For N = 2, the spin singlet sym- 
metry has been approximately restored with the Lipkin- 
Nogami approachi^^ Larger N have seldom been treated 
with the projection technique. ^'^ In view of the recent dis- 
cussion of spin effects in the transport spectra of quantum 
dots, information about the total spin is, however, neces- 
sary. Additionally, by restoring the symmetries correla- 
tions are introduced into the ground state wave function 
that are absent in a single UHF Slater determinant. This 
leads to a better estimate for the ground state energy. 
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In this paper, we apply a projection technique to the 
states obtained by UHF for estimating the ground state 
energy of a circular quantum dot with < 12 electrons, 
including a magnetic field. Starting from an UHF Slater 
determinant with broken rotational symmetry, a first es- 
timate for the ground state energy and the wavefunc- 
tion is obtained. Then, both the total spin and the an- 
gular momentum of the UHF variational wave function 
are introduced by projecting on the corresponding sub- 
spaces. We show that, after restoring all of the symme- 
tries, the energies and the wave functions are improved 
and show physical features which are not included in the 
UHF method. 

We discuss the efficiency of the projected HF method 
(PHF) by comparing our results with those of ED, CI, 
and QMC. We determine the ground state energies as 
a function of a magnetic field, and obtain the chemical 
potential that can be measured in transport experiments. 
Our main findings are: 

(i) By projecting the UHF wave functions on the to- 
tal angular momentum L and on the total spin 5*, the 
ground state energy is successively lowered. The correc- 
tion due to the spin projection is generally smaller than 
the one associated with the angular momentum, but still 
necessary for determining the correct ground state and 
its quantum numbers. 

(ii) The quantum numbers L and S are correctly repro- 
duced, if the strength of the interaction is not too large. 
Especially, for _B = 0, the first Hund's rule — namely 
that S is maximized for open shells — is recovered for 
iV < 12 electrons, except for N = 10, discussed below. 
Hund's rule has been claimed earlier to be violated on 
the basis of UHF results^. 

(m) By comparing the results with CI and QMC, we 
estimate a correlation energy, defined as the difference 
between PHF and "exact" energies, of about 2% of the 
ground state energy. 

{iv) With increasing interaction strength the correla- 
tion energy decreases. Nevertheless, for stronger inter- 
action, and larger TV, the PHF ground state tends to be 
spin polarized in contrast to more exact results. This 
is consistent with earlier conjectures, namely that UHF 
tends to overestimate the infiuence of the exchange.^^ 

[v) As a function of B, several crossovers between 
ground states with different total spins and angular mo- 
menta are found that are absent in UHF. These are asso- 
ciated with characteristic changes in the electron densi- 
ties. The onset of the singlet-triplet transition^^ occur- 
ring for dot filling factor v 2 and N even is recovered. 
Features that lead to an intrinsic spin blockade are pre- 
dicted. 

In the next Section, details of the UHF method are 
outlined. The consequences of the broken symmetries are 
described and the projection technique is discussed, with 
special emphasis on the total electron spin. In Sect. IIIII 
results for zero and non-zero magnetic field are presented 
and discussed. 



II. MODEL AND METHOD 
A. The model 

Consider N electrons in a two-dimensional (2D) quan- 
tum dot confined by an in-plane harmonic potential and 
subject to a perpendicular magnetic field B = Be^. The 
Hamiltonian is {h = c = 1) 



N ^ N 



(1) 



with 



[p + eA(r)] ^!!Mr^^g*,,Bs. (2) 
2m* 2 



r = (r, I?) the 2D polar coordinates, v{r) — e^/47reo£r'' 
the Coulomb interaction potential, B = rotA, m* effec- 
tive electron mass, luq confinement frequency, g* effective 
g-factor and /is the Bohr magneton. The z-component of 
the i-th spin is Szi — ±1/2, — e the electron charge and 
Eq (Sr) the vacuum (relative) dielectric constant. The 
single-particle term in ^ yields the Fock-Darwin^"' (FD) 
spectrum 

en,i,s^ =n{2n+ \l\ 1) + y ; + g*fiBBs, , (3) 

with eigenfunctions 0„,i(r)x^, where (x ) is the 
spinor corresponding to s^ = -1-1/2 (s^ = —1/2) and^ 



'?^'n,z(r) 




£o ^7r{n+\l\)l 

m _ ,.2 



(4) 



Here, n and I are principal and angular momentum quan- 
tum numbers, £q = (m*n)~2 the characteristic oscilla- 
tor length and ^i''(x) the generalized Laguerre polyno- 
mial. The cyclotron frequency ujc = eB/m*, and the 



effective confinement frequency fl — (w. 



?/4) 



1/2 



are 



introduced. 

At _B = 0, expressing energies in units ujq and lengths 
in units £q, the Hamiltonian ([T]) depends only on the 
dimensionless parameter 



(5) 



which represents the relative strength of the interaction. 

B. The unrestricted Hartree-Fock method 

In HF the Schrodinger equation for a given value of 
total Sz = Szi -I- ... -I- SzN is solved by using orbitals 



V'r(r) = <(r)x", 



l<i< Nry 



(6) 



3 



with a = + {a = —) denoting spin up (down) and Na is tained as the solutions of the coupled integro-differential 
the number of electrons with spin a ■ 1/2. They are ob- equations 



/io(r,al/2)+ / dr'p(r')w(r - r') 



C(r)-E 



drVr(r')C(r')«(r-r') 



^,"(r)=£rC(r), 



(7) 



where p{r) — p^(r) + p (r) is the HF density 



Na 



."(r)=^|V;r(r) 



(8) 



Here F^^^ are the Fock matrices, 



(14) 



i=l 



K 



For a given Sz, an initial guess for the orbitals "0" (i") 
with i < Na is made. Then, HF densities are evaluated 
and Eqs. ([7]) are solved to obtain updated orbitals. This 
is iterated until self-consistency is achieved. The many 
body wave function is a single Slater determinant, eigen- 
function of 5"^, 



1 



det{V'+(rj ),..., V^^(rj), 

V'r(rj),-.-,V';^_(rj)}, (9) 



A,i)=l 

and the two-body interaction matrix elements 
(^?7|w|i^A) 

dridr20* (ri)(/>*(r2)u(ri - r2)0^(ri)0A(r2) (15) 
can be evaluated analyticallyi^ The energy pO)) is then 

K 



that corresponds to a stationary point of the UHF en- E^' = ^ ^ ^ ai/2 f^i^M + -^"mI ^m'' • 
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(10) 



In order to numerically solve ^ we expand the orbitals 
in the FD basis (l)f,{r)x" (Eq. (g])) 



K 



l<f^<K (11) 



where C^j are complex coefficients. The truncation of 
the basis to K states is necessary in order to numerically 
implement the procedure. We have used the K = 75 
lowest FD states for each value of B. This led to fair 
convergence (see Sec. Ill E")) . Introducing the density ma- 
trices 



Q — ± fJ,,U—l 

We use spatially unrestricted initial conditions^i^iSl 
with a random distribution of initial C^^. This im- 
plies initial orbitals without circular symmetry, and leads 
to better energy estimates. However, symmetry broken 
Slater determinants are in general neither eigenfunctions 
of the total angular momentum L = li + . . . + In nor of 
(total spin S Si -I- . . . H- s^).^^ 
The most general UHF solution is a linear superposi- 
tion of eigenfunctions \'^^' {L, S)) of L and S 



oo N/2 

E 

L=-oo S>\S^ 



^ \^'^{L,S)). (17) 



Na 



connected to © by 



(12) 



For given N and Sz, many initial conditions are used. 
Correspondingly, several stationary points are found. 
They form a sequence (fc = 1,2,...) with ener- 



i=l 



gies Ef' < E2' < .... For a given 5*^, the process is 
iterated until the lowest E^' is found. The UHF ground 
state is defined as. 



K 



it is possible to show that equation ([7]) is equivalent to 
the coupled nonlinear Pople-Nesbet eigenvalue problem 



-EuHF = min |£;f" | . 



K 



^u^ui ^ fj,i 



(13) 



C. Spin and angular momentum projection 

In order to obtain states with specific L and S we act 
on the UHF Slater determinant with operators^^ Pl and 
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Pg' which project on L and S, respectively. They satisfy 
commutation rules [P|% Ft] = [Ps%H] = [Pl,H] 
Their simultaneous actii 
and S2, PiP|^|*«-) = 
energy is (Appendix E]) 



0. 

Their simultaneous action yields an eigenfunction of L 
and S^, PlP^^I^^')^ \^^'{L,S)). The corresponding 



The spin projector 



{^^'{L,S)\H\^'''{L,S)) 
(vl/S.(L,5)|vI/S.(L,^)) 



N/2 

n 



S'^-k{k + l) 
S{S + l)-k{k + l) 



(18) 



annihilates all the components of ([9]) with spin different 
from Its action is written as^i^ 



= Y,C,{S,S.,N)\T,^ 

9=0 



(19) where iV< = min{A^+, Af_} and 



(20) 



C,{S,S,,N) 



2S+1 



1 + N/2 + S 



s-s^ 

EM) 



q+S-S^-k 



k=0 



S S2 

k 



r 



s + s, 

S S 2, " k 



N/2 + S 
VN/2-q + k 



(21) 



are the Sanibel coefficientspi^'Ze The term \Tq) = \Tq 



(i)\ 



+ iTg^"'-*) is the sum of all 



(22) 



Slater determinants obtained by interchanging, without 
repetition, all the possible q spinor pairs with opposite 
spins in \^^'). By definition \Tq) = \^^'). 

For example, consider = 4, 5"^ = {N+ = iV_ = 2), 



1 



dct{Mj^X , ""2 X , "1 X ,"2X }• (23) 



This state is a linear superposition of all spin eigenstates 
with S <2. The spin projection selects a specific spin 



^s=o 




^") = 


1 

3 


\To) 


1 

~ 6 


|Ti>- 


1 

^3 


IT2), 


(24) 






^") = 


1 

2 


\To) 


1 

~ 2 








(25) 


pS,=0 




^") = 


1 

6 


\To) 


1 

+ 6 


\Ti)- 


1 

^6 


IT2), 


(26) 



where 
l^i) 



IT, 



[det{uj^X^,wJx ,'«iX^,'"2X } + 



+ det{u{x i'u^X^,uTx^,U2X } + 
+ det{u^'x~,W2'X"^,wrx"^,W2^X"} + 
+ det{u^x~,"2'x'^,"rx~,'"2'X^}] , 

^ =det{uj^X","Jx",^^rx^,^^2"X^} ■ 



Summing up equations ((24 



determinant l*?^), since J2s^S 



261) results in the original 



The projector on L is given by^ 

^27r 

Pl 



1 

2^./o 



d7 e-'^'>e'^'' , 



(27) 



where exp (iLj) acts on \Tq) rotating by 7 around the z 
axis all spatial parts of the orbitals 

<(r,79)^<(r,i? + 7). 

We denote this by |r5(7)). Using ^ and ^ we get 

27r 



1 ^ 

|*^iL,5)) = ^C,(5,5„iV). 



g=0 
27r 



d7 e-^^^|T,(7)) 



(28) 



The projected state ([^5]) is a sum of many Slater deter- 
minants (Appendix [X)) . This indicates that correlation 
has been introduced by the projection. 

The main computational effort is due to the evaluation 
of two-body matrix elements in p^ . Projecting an A^- 
particle UHF state with 5*2 to a state with total spin S 
requires to evaluate n{Sz, N) — X]^o terms, 

N 

^N/2 + S, 



(29) 



For A^ even (odd), the worst case is 6*2 = {Sz = 1/2). 

For the angular momentum projection we use a fast 
Fourier transform (FFT) and partition the integration 
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interval [0, 2tt\ in n{L) points; n{L) is determined by the 
angular momentum range \L\ < L^^x for which good 
convergence (relative error < 10~^) of the PHF ener- 
gies is required. We have checked that for L^ax = 20 
n(L) — 256 is needed. Using FFT, all energy values for 
given S and \L\ < Lmax are simultaneously available, 
which considerably accelerates the calculation with re- 
spect to performing distinct computations for each value 
of L. The total number of two-body matrix elements is 
n-tot = n{S^,N)n{L). 



N_ ntot n_ 

2 512 

4 1536 19774^" 

6 5120 661300^^ 

8 17920 

10 64512 

12 236544 

14 878592 

16 3294720 



TABLE I: Numbers of two-body matrix elements used for 
PHF, ntot required to evaluate E^'^^{L, S) for even A'' < 16, 
Sz = S < N/2, and \L\ < 20 (see text). Last column: some 
numbers used in other methods. 

Table U shows ntot for the case of even iV < 16 and 
Lmax = 20 in the worst case Sz = 0. This is the Lmax 
value used in the paper. Although ntot quickly increases 
as a function of N , especially because p9|) grows expo- 
nentially for large iV, it still compares favorably with 
respect to exact methods. For example, previously re- 
ported ED calculations'^'' used a basis of 19774 Slater 
determinants for = 4, ^2 = 0, 5' = 2, and L = 14. 
CI calculations^^ for iV = 6, 5'^ = 0, 5* = 0, L = need 
661300 configurational state functions (linear superposi- 
tion of Slater determinants). 

D. Determining the PHF ground state 

To determine the ground state, it is generally not suf- 
ficient to project only the UHF ground state. If several 
UHF solutions (Sec. IIIB[) are almost degenerate, all of 
the l^'f^) have to be projected. The PHF ground state 
is then defined by 

i;PHF= min {Ef^{L,S)] . (30) 

One can show that projecting arbitrary UHF Slater de- 
terminants on L and S always leads to energies that are 
not lower than the exact ground state energy, thus satis- 
fying the variational principle. 

As an example, we consider = 4 for i? = 0, with con- 
finement energy wq = 0.741 meV. We assume the stan- 
dard GaAs parameters, m* = 0.067 rrie and = 12.4. 



The confinement corresponds to A = 4. For each Sz, 
we have used more than 500 UHF initial conditions. We 
found two solutions with Sz = 0, and one for 5*2 = 1 
and Sz — '2- The corresponding energies Ef'' are given 
in Tab. |TT1 



Sz 




B^(L,S) 


L S 





19.612 


19.356 









19.404 


1 1 




19.641 


19.331* 


1 






19.515 


2 


1 


19.608 


19.342 


1 






19.394 


1 1 


2 


19.581* 


19.516 


2 2 



TABLE II: PHF for A = 4 for a GaAs quantum dot with 
B = 0, Wo = 0.741 meV, m* = 0.067 me, £r = 12.4, A = 
4; 2nd column: UHF energy Ef' (units uq); 3rd column: 
PHF energy Ef~ (L, S) (units uq); 4th column: PHF quantum 
numbers. Ground states are indicated by *. 

The UHF ground state corresponds to S'^ = 2. Ap- 
plying the projection to each UHF state, we obtain 
Ef'{L, S) with different quantum numbers L and S: the 
lowest two are given. The PHF ground state corresponds 
to L = and S = 1. The latter is obtained by projecting 
the energetically higher UHF minimum with Sz = 0- 



E. Some comments about errors 

The major systematic error of the UHF approximation 
is the neglect of correlations. By projecting the slater 
determinant on fixed angular momentum and spin PHF 
attempts to correct for these effects. A second systematic 
effect is due to the uncertainty if the self consistent HF 
procedure has converged towards the absolute minimum 
of the energy. 

In determining the UHF ground state energies, we have 
checked that the convergence with respect to the size of 
the basis set is better than 10~®. 

For getting insight into the above systematic effects 
one can start from wave functions with the same L, S 
but originating from UHF states with different Sz ■ They 
should be degenerated at i? = 0. In the example of 
Tab. El these are the pairs |*?(1,1)), |*}(1,1)) and 
1*^(0,1)), |^'}(0,1)). Their energetic differences are 
0.010 ojQ and 0.011 ojo, respectively. This corresponds 
to a relative uncertainty of 5 • 10~*. Similar estimates for 
the "degeneracy error" is obtained from data for different 
N and A. We attribute the degeneracy error mainly to 
UHF: different UHF states in different Sz sectors approx- 
imate the true states with different precision. Therefore, 
their projection on the same L, 5* sector does not yield 
exactly degenerate states. 

By comparing our results with other works (see below). 
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the PHF ground state energies for N > 3 remain about 
2% higher than those obtained with ED and QMC. This 
can be attributed to correlations beyond those introduced 
by the projection. This also is the limiting factor for the 
ground state quantum numbers in the regime N > 6 and 
A > 4 where too high polarization are obtained. 

When several PHF energies Ef'{L,S) are almost de- 
generate, one can improve further the ground state: 
linear superposition of the almost degenerate states 
|\E'f^(L, S)) may result in further lowering of the energy. 
Here, we have not systematically investigated this effect. 

III. RESULTS 
A. Zero magnetic field 

1. Ground state energies 



N 


A 


i?PHF 


L S 


Eci 


-Bdmc 


L S 


2 


1.89 


3.817 







3.649 







2 


3.885 





3.7295 









4 


4.983 





4.8502 







3 


1.89 


8.154 


1 1/2 




7.978 


1 1/2 




2 


8.337 


11/2 


8.1671 




11/2 




4 


11.131 


3/2 


11.043 




1 1/2 


4 


1.89 


13.554 


1 




13.266 


1 




2 


13.899 


1 


13.626 




1 




4 


19.330 


1 


19.035 




1 


5 


1.89 


20.264 


1 1/2 




19.764 


1 1/2 




2 


20.811 


1 1/2 


20.33 




1 1/2 




4 


29.501 


1 1/2 


28.94 




1 1/2 


6 


1.89 


27.905 







27.143 







2 


28.703 





27.98 









4 


41.187 


3 


40.45 







7 


1.89 
2 


36.627 
37.698 


2 1/2 
2 1/2 




35.836 


2 1/2 




4 


54.497 


5/2 


(54.68) 


53.726 


(2)2 (1/2)1/2 


8 


1.89 


46.260 


1 




45.321 


1 




2 


47.659 


1 


47.14 


46.679 


1 




4 


69.479 


4 


70.48 




1 


9 


1.89 


56.853 


3/2 




55.643 


3/2 


10 


1.89 


68.245 







66.8785 


2 1 






(68.283) 


2 1 




(66.8789) 





11 


1.89 


80.444 


1/2 




78.835 


1/2 


12 


1.89 


93.661 







91.556 






TABLE III: Ground state energies from PHF for TV < 12 and 
A = 1.89, 2, 4 with corresponding L, S (m* = 0.067 me, and 
= 12.4) together with results from CI^^, and DMC (Ref.|39| 
for A — 1.89, Ref. [3^ for A > 2). All energies are in units loq- 

Table IIIII summarizes our results for the ground state 



energies at B = 0, for N < 12 and A = 1.89, 2, 4. 
Results obtained with Diffusion Monte Carlo (DMC)^i2a 
and CI^ are included. 

For A < 2, angular momenta and total spins of the 
ground states obtained by PHF agree with DMC and CI. 
The total spin fulfills Hund's first rule: a singlet state 
for the filled shells {N = 2, 6, 12), a triplet for = 4, 8 
and S = 3/2 for N = 9. Only for A^ = 10, Hund's rule 
is not fulfilled since we find S = instead of 5 = 1. 
However, here the degeneracy error is 0.064 wq, larger 
than the energy distance AE = 0.038 wq between the 
ground and the first excited state. Also DMCS2, predicts 
an extremely small energy gap between the singlet and 
the triplet, though it yields an 5* = 1 ground state. 

Increasing the interaction strength (A — 4) PHF still 
produces energies consistent with CI and DMC. However, 
for A^ — 3, 6, 7, 8 incorrect quantum numbers are pre- 
dicted with a tendency towards polarization. Whenever 
polarization occurs the ground states have low angular 
momenta in PHF. 

For A > 4, preliminary results indicate deviations of 
PHF with respect to CI, DMC. They are reminiscent 
of the tendency of HF to predict spin polarized ground 
states due to overestimating the exchange as compared 
to correlations. 

The relative deviation (5 — {EpHp— Edmc) / Enuc for— 
A = 1.89 and 2 < A^ < 12 is shown in Fig. [11 (5 is largest 
for the closed shells A^ = 2,6,12. Except for N = 2, 
(5 w 2%. The inset shows S for N = 2 (squares, with 
L, S" = 0, 0), A^ = 4 (dots, with L,S ^0, 1), and A^ = 6 
(triangles, with L,S = 0,0) within 1.89 < A < 8. A 
decrease with A according to a power law is observed, 
\-l3{N) _ gy numerically fitting the data, one finds /3(2) — 
0.57, /3(4) = 0.44 and /3(6) = 0.45. 




FIG. 1: Deviations 5 between PHF and DMO^ for 2 < A < 
12, with A = 1.89 (Tab. |ITI1). Inset: double logarithmic plot 
of (5(A) from PHF and QMGS^ (A = 1.89), QMC''** and GI^ 
(A > 2) for A = 2 (squares), A = 4 (dots) A = 6 (triangles). 
Lines: best fits to data. 

Table ITVl and Fig. [5] illustrate the effect of angular mo- 
mentum projection alone followed by spin projection for 
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Ef^iL) 
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i^DMC 


48.150 
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47.659° 
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48.031^ 
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46.875 
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47.971 
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48.028^ 
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48.025 
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47.910^ 
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48.131 
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47.796 
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47.742'' 
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46.756 
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46.917 
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48.022 
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47.997° 
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48.243 
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47.896 
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47.881^ 
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48.335 
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48.129 
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48.126 
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47.404 



TABLE IV; Comparison of the lowest energies Ef' obtained 
from UHF, followed by projection on angular momentum, 
Ef'{L), and total spin, Ef'{L,S) for iV = 8, A = 2 (m* = 
0.067 me, Er = 12.4). Last column: energies from DMCj^^ 
The ground state has L, S = 0, 1. Superscripts 0'i'2,3,4 (jgjjQ^g 
"degenerated" energies with the same quantum numbers L, S 
but originating from different Sz- All energies are in units loq. 



El = 48.150 




E^(Q,l) = 47.659 



UHF L-projection PHF 

FIG. 2: Influence of the projection procedure on the energy 
levels (unit uoq) for A'' = 8 and A = 2. Only the two lowest 
UHF states (left) directly involved in the determination of the 
PHF ground state (right) are shown. 

state, in contrast to the correct result. 

The degeneracy error for this case is approximately 
0.06 (some example of almost degenerate states are 
included in Tab. lIV[) . The distance between ground state 
and the first excited state is w 0.08 wq. This suggests 
that the ground state for = 8 has L,S = Q, 1, consis- 
tent with DMC. Even the quantum numbers of the first 
three excited states turn out to be reproduced correctly 
while the 4th and the 5th appear to be interchanged. 



A^ = 8 starting with UHF states with = 0, . . . , 3. The 
energy Ef" (L) projected on angular momentum is 



{^f'\H\^f'iL)) 



(31) 



where |4'f"(L)) Pil^f"). Only the lowest energies are 
included in the table. 

The typical energy gain obtained by angular momen- 
tum projection is about 0.25 o^q. The spin projection in- 
duces corrections of the same order of magnitude, which 
can even change the sequence of energies (Fig. [5]). From 
the UHF state with 3^ = and = 48.150^0, which 
is not the UHF ground state, projection on L = 
yields El{L = 0) = 47.842 cjq- After projection on the 
total spin we obtain the energy of the ground state, 
Ei{L,S = 0,1) = 47.659 and an excited state at 
Ef{L, 5 = 0, 2) = 48.031 wq. On the other hand, the en- 
ergetically lowest UHF minimum El = 48.131^0 turns 
out to yield the first excited PHF state at E\{L,S = 
2,1) = 47.742 Wo- 

Thus, PHF not only introduces a lowering of the en- 
ergies but can also restore the correct ordering of energy 
levels. This can be seen from the last column of Tab. lIVi 
which contains the results obtained by DMG^^. Restor- 
ing the spin plays a crucial role in obtaining all correct 
quantum numbers for the ground state including Hund's 
rulci^ For example with angular momentum projection 
alone, one would have predicted L = 2 for the ground 



2. Ground state densities 



For the spin-resolved densities 



Pphf('^) 



i^'' (L, g)| S{r - r.)^.,„c..i/2l^^- {L, S)) 



we first consider A^ = 3 and A^ = 4 (Figs. E] and S]) for 
intermediate (A = 2) and strong (A — 8) interaction. 
Increasing the interaction strength leads to a shift of the 
maximum of the densities towards higher r, consistent 
with earlier findings by ED.^-i^ This is clearly observed 
in the spin up density for A^ = 3 (Fig. [3]). For A^ = 4, the 
ground state (L = 0, S* = 1, 5'^ = 0) densities Ppjjp(r) = 
Pphf('') (I^ig- SI agree very well with ED^^ for large r. 
Generally, deviations occur near r « 0. 

Figure [5]^a) shows the total electron density for N — 5, 
L,S — 1,1/2, for A — 0.5,2,10. For weak interaction, 
A — 0.5 (solid line), we find good agreement with CI — 
(squares). For A = 2 (dashed) small deviations near r = 
are found. Figure (Hb) indicates that the spin-down 
density is responsible for the small deviation from the 
exact result around r = for A = 2. 



B. Finite magnetic field 

In this section, we show results for A^ = 5, 6, 7 in the 
presence of a magnetic field, B < 2.4 T, corresponding 
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FIG. 3: Spin resolved densities Pphf(7') (thick line: a = +; 
thin line: a = —) for a GaAs quantum dot with A'' = 3, 
L = 1, S = 1/2, and Sz = 1/2 for interaction strengths \ — 2 
(solid), A = 8 (dashed). Density unit: it~^£q'^. Data from 
ED.>^ squares A = 2, circles A = 8. 



0.5 




12 3 4 

r 

FIG. 4: Spin resolved densities Pphf('') ~ Pphf('') f^'' ^ GaAs 
quantum dot with A'^ = 4, L = 1, S = 1, and Sz = for A = 2 
(solid line), A = 8 (dashed line). Density unit: tt'^Iq"^ . Data 
from ED:2^ squares A = 2, circles A = 8. 



to a dot filling factor v ^ 2 {N < 4 has been discussed 
in Ref . 'ts') . We assume here a confinement t^o = 6 meV 
(corresponding to A = 1.45) and g* = —0.44. For B > 0, 
due to the Zeeman term, the PHF ground state always 
has 5*2 = S. Therefore, we do not specify Sz in the 
following. 

We start with = 5 (Fig. ED and A^ = 6 (Fig. E]). We 
show the UHF ground state energy -Euhf (solid line), 
the energy obtained from angular momentum projection 
(dashed, Eq. (pij) '). and the PHF energy (dashed-dotted, 
Eq. (PO)) ). The highest energy gain is here due to the 
angular momentum projection. Spin projection leads to 
a further decrease of the ground state energy. Obviously, 
UHF and PHF results behave completely differently with 
B. 

For instance, for A^ = 5 (Fig. [5]) the UHF ground state 
shows crossovers 5^ = 1/2 — > 3/2 at i? sa 0.9 T and Sz = 




r 4 r 4 



FIG. 5: Densities for A^ = 5 L = 1, S = 1/2, and Sz = 
1/2. (a) Total density pphf(j") = Ppupir) + Ppjjp(r) (units 
TT~'^£o^) for A = 0.5 (solid), A = 2 (dashed), A = 10 (dashed- 
dotted). Squares, circles and triangles: data from CI.— (b) 
Spin resolved densities pPhf('') (solid: a = +, dashed: a = —, 
units tt'^Eq^) for A = 2. Squares, circles: data from CI. (c) 
Same as (b) but A — 10. 
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0.6 1.2 1.8 2.4 

B 

FIG. 6: Ground state energy E (units loo) as a function of 
magnetic field B (units T) for A^ = 5. Solid: UHF; dashed: 
angular momentum projection; dashed-dotted: PHF. Here, 
and in the following figures, m* — 0.067 me, Sr ~ 12.4, g* = 
—0.44, and uiq = 6 meV. 



3/2 ^ 1/2 at B w 1.5 T. In contrast, the PHF energy 
has total spin S = 1/2 = Sz in the entire magnetic field 
region. The state with Sz = 3/2, not compatible with the 
total spin S = 1/2, is certainly an artifact of UHF. The 
crossover L, 5" = 1, 1/2 ^ L, 5 = 4, 1/2 with increasing 
magnetic field at B = 1.4 T agrees quantitatively with 
the earlier results obtained by ED.^^ 

For Af = 6 (Fig. U}, UHF (sohd Hue) displays no Sz 
transitions. When rotational symmetry is restored, two 
crossovers, (a) L = ^ — 3 and (b) L = — 3 ^ — 6 ap- 
pear. Performing the spin projection, singlet states cor- 
responding to L = and L = — 6 are found, and 5 = 1 
for L = —3 is obtained. Singlets have the largest energy 
gain, leading to a shift of the features found with angu- 
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FIG. 7: Ground state energy E (units uo) as a function of 
magnetic field B (units T) for TV = 6. Solid: UHF; dashed: 
angular momentum projection; dashed-dotted: PHF. 



32.15 



31.60 




L=-2 S=l/2 



FIG. 8: Ground state energy E (units luo) as a function of 
magnetic field B (units T) for TV = 7. Solid: UHF; dashed: 
angular momentum projection; dashed-dotted: PHF. 



lar momentum projection (Fig. [7|). Also here, the PHF 
quantum numbers agree with the earlier results obtained 
by ED/^ including the magnitudes of the crossover fields 
at B w 1 T and B « 1.8 T respectively. 

The singlet-triplet crossover occurring for TV = 6 at 
B « 1.8 T, corresponds to a filling factor v k, 2, and is 
a peculiar feature which is confirmed by several experi- 
mental and theoretical studies J^ii^i^iSS Also for TV = 8 
preliminary data indicate such a crossover near v = 2. 
These crossovers are completely absent in UHF (Fig.[71). 

Most interesting is TV = 7 (Fig. [8]): near B w 2.2 T 
the ground state has S — 3/2. This can only be obtained 
including the spin projection and leads eventually to a 
spin blockade in the transport (see below). 

In Fig. [9] we show the scheme of the ground state quan- 
tum numbers for 4 < TV < 7, as obtained by PHF. They 



31.79 



per- 



qualitatively agree with previous calculations 
formed for TV < 6. In the region of B, where the ground 
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FIG. 9: Scheme of the quantum numbers L, S of the PHF 
ground state as a function of the magnetic field B (units T) 
for 4 < AT < 7. 



state of TV = 7 has S — Z/2 the state with TV = 6 is a 
singlet. Since AS* > 1/2 between the two ground states 
a spin blockade in the 6 7 transition can be expected 
near the edge oiv — 2 for N = 7 electrons, for B « 2.3 T. 
We note in passing that the lowest excited states for 
TV = 7, with L = -5,S = 1/2 and L = -9,S = 1/2, 
are at most w 0.07 meV (~ 0.8 mK) higher in energy. 
Therefore, it may be difficult to experimentally observe 
this blockade. 

The chemical potential traces ^n{B) ~ £'phf(^, B) — 
£'phf(A^ — 1, B) obtained by PHF when varying B are ex- 
perimentally accessible via Coulomb blockade. Figure [TUl 




FIG. 10: Chemical potentials ^^[B), ^q[B), fJ.7{B) (units ujq) 
as a function of B (unit T). Arrows: edge of filling factor 
= 2 for TV = 5 (bottom panel), TV = 6 (center), TV = 7 
(top). Red line: region of intrinsic spin blockade (see text). 

shows /i5(i3), neiB) and ij.t{B). Arrows indicate the on- 
set oi V = 2 for the configuration with TV = 5 (bottom 
panel), TV = 6 (center), N — 7 (top). The chemical po- 
tentials exhibit features related to the above discussed 
crossovers between ground states. At the onset oi v — 2, 
the chemical potentials exhibit a cusp. For even TV, this 
corresponds to the above mentioned singlet-triplet tran- 
sition. --4 Generally, the chemical potentials show kinks 
when quantum numbers of the ground states change 
(Figs, [i and E 
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IV. CONCLUSION 



In order to obtain (|T8l) , we have to evaluate the overlaps 



We have described a systematic procedure to overcome 
some of the limitations of UHF approach. Using angu- 
lar momentum and total spin projections, we have intro- 
duced correlations that provide lower estimates for the 
ground state energies, besides determining the spin and 
the angular momentum. Several sources of errors have 
been discussed. In particular, a degeneracy error has 
been found to be useful for deciding whether or not the 
estimate for the ground state is plausible. 

The procedure yields results consistent with earlier 
findings for interaction strengths A < 2 which corre- 
sponds to experimentally relevant confinement energies 
Wo > 3 meV for = 12.4,^ 

For B — Q and A < 2, we have confirmed Hund's first 
rule for the dot total spin, except for N = 10. In this 
case, the ground state is ambiguous, since the energy gap 
between ground and first excited state is smaller than 
the degeneracy error, consistent with other results. For 
stronger interaction, A > 4, deviations from Hund's rules 
are obtained, accompanied by the well-known exchange 
induced tendency of HF-based methods to favor ground 
states with higher spins and zero angular momenta. 

We have shown that PHF predicts correctly the fea- 
tures of the ground state energy as a function of B. 
We have found a spin blockade in the transport between 
N = 6 and N = 7, occurring at a filling factor « 2. 

Given the slower increase in computational effort with 
particle number described in Sec. Ill CI (Tab. [l|, as com- 
pared to other methods, we hope by parallelization of our 
code to obtain in the future results for higher number of 
particles {N > 20), varying B, for interaction strengths 
relevant to quantum dot experiments, A < 2. 

That the densities are correctly reproduced suggests 
that tunneling rates between the quantum dot and at- 
tached leads needed for electron transport can be rea- 
sonably well estimated when using PHF wave functions. 
This might be useful for providing quantitative results 
for predicting the heights of the Coulomb blockade peaks 
as a function of _Bjiii^ii^ 

This work has been supported by the Italian MIUR via 
PRIN05, by the European Union via MRTN-CT-2003- 
504574 contract and by the SFB 508 "Quantenmateri- 
alien" of the Universitat Hamburg. 



APPENDIX A: SOME DETAILS OF THE 
IMPLEMENTATION 



We provide some technical details about the implemen- 
tation of the projection technique outlined in Sec. Ill CI 



^^C,(5,5.,iV)(To|T«(7)), (Al) 

q=0 i=l 

and the Hamiltonian matrix elements 



(^■^^ \Ho + (L, S)) = — d7 e-'^-^ • 

g 5., N){T,\H, + y|TW(^)) , (A2) 

9=0 1=1 

with 



N 



1 ^ 



(A3) 



In the following, we avoid explicit reference to electron 
coordinates unless when strictly necessary. 

Previously, analytic expressions for (|18p have been re- 
ported^^ for N — 2. For larger N, however, the num- 
ber and the complexity of the above matrix elements in- 
creases dramatically. This eventually requires numerical 
treatment. 

For evaluating |T,) = Y.Zi \Tg''^): 

we need to generate 
all the Uq swaps of q opposite spin pairs in |To). The lat- 
ter correspond to a special class of permutations, acting 
on the fc-th component of the generalized vector 



- = (x+, 

with the correspondence k ~f n^. 
with 



■,X-) (A4) 
One then has 



WAfCr } , (A5) 



(A6) 



All the permutations are pre-tabulated at the beginning 
of the calculation. Further calculations are performed 
by means of well known theorems^ for many body wave 
functions. The overlap term is 

(ro|rW(7))=det{d('-^)(7)} (A7) 

where d^''''-' is the overlap matrix 

4f (7) = Kkp(7)>K>4,.)). (AS) 

Here, 1 < k,p < N and Wp{-f) is a shorthand nota- 
tion for the rotated spatial part Wp{r,-d + 7). The term 
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(crfc|(T^(,,i)) reduces to a Kronecker delta. For the evalua- 
tion of {wk\wp{'y)) we use the FD basis, which is particu- 
larly convenient to describe rotations. The spatial parts 
transform as 

K 

^ + 7) = E C;re*'-^0^(r, 0) (A9) 
with the angular momentum of the /i-th FD state, 



for p < iV_(- and ap 



for N+ + 1 < p < N. 



Therefore, 



7 



(AlO) 



In the interaction part 



1 ^ 



N 

2 / ^ / ^ feife2PlP2 



felfe2|PlP2 



(7) 



(A13) 



^fcife! |piP2 '•^•^ represents the second order cofactor of the 
matrix d(«'*)(7), and 



which can easily be evaluated and stored. 

The single particle term in the Hamiltonian is 



fclfc2PlP2 ^ ' '' 



N 

(To|i7o|Tf (7)> - E 4'fe?(^)^i|P^(^) (All) (u'fe,^«,J«|z«,,(7)z«p2(7)>(<TfeJa^<,,o)(a,2k4,.)) 



1 



where -DI,?'*'' (7) is the fc,p first order cofactor of c?^'^'*^ and 



■k\p 

*&p(7) = (wfc,crfe|/lo|Wp(7),CT^(' 



(A12) 



Terms H^%'^{"f) and VfcXp (7) can be straightfor- 
wardly evaluated as for (jA10|) . 



^ See e.g. S. M. Reimann, and M. Manninen, Rev. Mod. 
Phys. 74, 1 283 (2002); M. Manninen and S. M. Reimann, 
|cond-mat/0 703292 , and references therein. 

^ I. Romanovsky, C. Yannouleas, and U. Landman, Phys. 
Rev. Lett. 93, 230405 (2004); L. O. Baksmaty, C. Yan- 
nouleas, and U. Landman, Phys. Rev. A 75, 023620 (2007). 

^ A. Bohr and B. R. Mottelson, "Nuclear structure", Lon- 
don, W. A. Benjamin (1975). 

* W. de Heer, Rev. Mod. Phys. 65, 611 (1993). 

^ L. P. Kouwenhoven, C. M. Marcus, P. L. McEuen, S. 
Tarucha, R. M. Westervelt, and N. S. Wingreen, Proceed- 
ings of the NATO Advanced Study Institute on Mesoscopic 
Electron Transport, L. L. Sohn, L. P. Kouwenhoven and 
G. Schon editors (Kluwer Series E345) p. 105 (1997). 

® L. P. Kouwenhoven, T. H. Oosterkamp, M. W. S. Danoe- 
sastro, M. Eto, D. G. Austing, T. Honda, and S. Tarucha, 
Science 278, 1788 (1997). 

S. Tarucha, D. G. Austing, T. Honda, R. J. van der Hage, 
and L. P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996); 
S. Sasaki, D. G. Austing, and S. Tarucha, Physica B 256, 
157 (1998). 

* R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, cond-mat/0610433 

® D. Weinmann, W. Hausler, and B. Kramer, Phys. Rev. 

Lett. 74, 984 (1995). 
° A. K. Huettel, H. Qin, A. W. HoUeitner, R. H. Blick, K. 

Neumaier, D. Weinmann, K. Eberl, and J. P. Kotthaus, 

Europhys. Lett. 62, 712 (2003). 
^ M. C. Rogge, F. Cavaliere, M. Sassetti, R. J. Hang, and 

B. Kramer, New J. Phys. 8, 298 (2006). 
^ M. Ciorga, M. Pioro-Ladriere, P. Zawadzki, P. Hawrylak, 



and A. S. Sachrajda, Appl. Phys. Lett. 80, 2177 (2002). 
F. Cavaliere, A. Braggio, J. T. Stockburger, M. Sassetti, 
and B. Kramer, Phys. Rev. Lett. 93, 036803 (2004); F. 
Gavaliere, A. Braggio, M. Sassetti, and B. Kramer, Phys. 

Rev. B 70, 125323 (2004). _ 

B. Muralidharan and S. Datta, cond-mat/0702i61 
M. Ciorga, A. Wensauer, M. Pioro-Ladriere, M. Korkusin- 
ski, J. Kyriakidis, A. S. Sachrajda, and P. Hawrylak, Phys. 
Rev. Lett. 88, 256804 (2002). 

M. Ciorga, A. S. Sachrajda, P. Hawrylak, C. Gould, P. 
Zawadzki, S. Jullian, Y. Feng, and Z. Wasilewski, Phys. 
Rev. B 61, R16315 (2000). 

M. C. Rogge, C. Fuhner, U. F. Keyser, and R. J. Hang, 
Appl. Phys. Lett. 85, 606 (2004). 

D. Kupidura, M. C. Rogge, M. Reinwald, W. Wegscheider, 
and R. J. Hang, Phys. Rev. Lett. 96, 046802 (2006). 
^3 M. Stopa, W. G. van der Wiel, S. De Franceschi, S. 
Tarucha, and L. P. Kouwenhoven, Phys. Rev. Lett. 91, 
046601 (2003). 

^° M. C. Rogge, C. Fuhner, and R. J. Hang, Phys. Rev. Lett. 

97, 176801 (2006). 
21 M. Keller, U. Wilhelm, J. Schmid, J. Weis, K. v. Klitzing, 

and K. Eberl, Phys. Rev. B 64, 033302 (2001). 

M. Dineykhan and R. G. Nazmitdinov, Phys. Rev. B 55, 

13707 (1997). 

D. Pfannkuche, V. Gudmundsson, and P. A. Maksym, 
Phys. Rev. B 47, 2244 (1993). 
'^'^ M. Wagner, U. Merkt, and A. V. Chaplik, Phys. Rev. B 
45, 1951 (1992). 

S. A. Mikhailov and N. A. Savostianova, Phys. Rev. B 66, 
033307 (2002). 



12 



^'^ S. A. Mikhailov, Phys. Rev. B 66, 153313 (2002). 

P. Hawrylak, D. Pfannkuche, Phys. Rev. Lett. 70 485 
(1993). 

P. Hawrylak, Phys. Rev. Lett. 71 3347 (1993). 
2^ A. Wojs and P. Hawrylak, Phys. Rev. B 53, 10841 (1996). 
^° M. B. Tavernier, E. Anisimovas, F. M. Peeters, B. Szafran, 

J. Adamowski, and S. Bednarek, Phys. Rev. B 68, 205305 

(2003). 

H. Imamura, H. Aoki, and P. A. Maksym, Phys. Rev. B 

57, R4257 (1998). 

J. Kyriakidis and C. J. Stevenson, cond-mat/0608044 
3^ Y. Nishi, P. A. Maksym, D. G. Austing, T. Hatano, L. P. 
Kouwenhoven, H. Aoki, and S. Tarucha, Phys. Rev. B 74, 
033306 (2006). 

A. Wensauer, M. Korkusinski, and P. Hawrylak, Solid 
State Commun. 130, 115 (2004). 

M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, J. 

Chem. Phys. 124, 124102 (2006). 
^® K. Varga, P. Navratil, J. Usukura, and Y. Suzuki, Phys. 

Rev. B 63, 205308 (2001). 
3'' A. Ghosal, A. D. Giiglii, C. J. Umrigar, D. UUmo, and H. 

U. Baranger, Nature Physics 2, 336 (2006). 
^® A. Ghosal, A. D. Giiglii, C. J. Umrigar, D. UUmo, and H. 

U. Baranger, arXiv:cond-mat/0702666 

F. Pederiva, C. J. Umrigar, and E. Lipparini, Phys. Rev. 

B 62, 8120 (2000); F. Pederiva, C. J. Umrigar,' and E. 

Lipparini, Phys. Rev. B 68, 089901(E) (2003). 

L. CoUetti, F. Pederiva, E. Lipparini, and C. J. Umrigar, 

Eur. Phys. J. B 27, 385 (2002). 
'^^ R. Egger, W. Hausler, C. H. Mak, and H. Grabert, Phys. 

Rev. Lett. 82, 3320 (1999); R. Egger, W. Hausler, C. H. 

Mak, and H. Grabert, Phys. Rev. Lett. 83, 462 (1999). 
''^ A. Harju, S. Siljamaki, and R. M. Nieminen, Phys. Rev. B 

65, 075309 (2002). 

F. Bolton, Phys. Rev. B 54, 4780 (1996). 
E. Rasanen, H. Saarikoski, A. Harju, M. Ciorga, and A 
Sachrajda, unpublished. 
^5 P. A. Maksym, Phys. Rev. B 53, 10871 (1996). 

W. Y. Ruan, Y. Y. Liu, C. G. Bao, and Z. Q. Zhang, Phys. 
Rev. B 51, 7942 (1995). 

A. V. Filinov, M. Bonitz, and Yu. E. Lozovik, Phys. Rev. 
Lett. 86, 3851 (2001). 

C. G. Bao, W. Y. Ruan, and Y. Y. Liu, Phys. Rev. B 53, 
10820 (1996). 

*® C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 5325 

(1999) ; C. Yannouleas and U. Landman, Phys. Rev. Lett. 
85, 2220 (2000). 

C. Yannouleas and U. Landman, Phys. Rev. Lett. 85, 1726 

(2000) . 

^'^ B. Reusch and H. Grabert, Phys. Rev. B 68, 045309 
(2003). 

^'^ B. Reusch, W. Hausler, and H. Grabert, Phys. Rev. B 63, 
113313 (2001). 



A. Emperador, E. Lipparini, and LI. Serra, Phys. Rev. B 
73, 235341 (2006). 

A. Wensauer, M. Korkusinski, and P. Hawrylak, Phys. 
Rev. B 67, 035325 (2003). 

P. Hawrylak, C. Gould, A. Sachrajda, Y. Feng, and Z. 
Wasilewski, Phys. Rev. B 59, 2801 (1999). 
M. Koskinen, M. Manninen, and S. M. Reimann, Phys. 
Rev. Lett. 79, 1389 (1997). 

K. Hirose, and N. S. Wingreen, Phys. Rev. B 59, 4604 
(1999). 

M. Gattobigio, P. Capuzzi, M. Polini, R. Asgari, and M. 
P. Tosi, Phys. Rev. B 72, 045306 (2005). 
^® A. Harju, H. Saarikoski, and E. Rasanen, Phys. Rev. Lett. 
96, 126805 (2006). 

E. Rasanen, J. Konemann, R. J. Haug, M. J. Puska, and 

R. M. Nieminen, Phys. Rev. B 70, 115308 (2004). 

H. Saarikoski, E. Rasanen, S. Siljamaki, A. Harju, M. J. 

Puska, and R. M. Nieminen, Phys. Rev. B 67, 205327 

(2003). 

N. Helbig, S. Kurth, S. Pittalis, E. Rasanen, and E. K. U. 
Gross, arXiv:cond-mat/0605599 

A. Szabo and N. S. Ostlund, "Modern Quantum Chem- 
istry" Dover (New York) (1999). 

P. -O. Lowdin, Phys. Rev. 97, 1509 (1955). 
P. Ring and P. Schuck, "The Nuclear Many-Body Prob- 
lem" Springer- Verlag (New York) (1980). 
P. -O. Lowdin, Phys. Rev. 97, 1474 (1955). 
C. Yannouleas and U. Landman, Phys. Rev. B 68, 035325 

(2003) . 

^® C. Yannouleas and U. Landman, Phys. Rev. B 69, 113306 

(2004) . 

Y. Li, C. Yannouleas, and U. Landman, Phys. Rev. B 73, 
075301 (2006). 

^° H. -M. Miiller and S. E. Koonin, Phys. Rev. B 54, 14532 
(1996). 

'^^ C. Yannouleas and U. Landman, J. Phys.: Condens. Mat- 
ter 14, L591 (2002). 

LI. Serra, R. G. Nazmitdinov, and A. Puente, Phys. Rev 
B 68, 035341 (2003). 

U. D. Giovannini, F. Cavaliere, R. Cenni, M. Sassetti, and 

B. Kramer, New J. Phys. 9, 93 (2007). 

V. Fock, Z. Phys. 47, 446 (1928); G. C. Darwin, Proc. 

Cambridge Philos. Soc. 27, 86 (1930). 
^5 V. H. Smith, J. Chem. Phys. 41, 277 (1964). 
"^^ M. B. Ruiz, J. Math. Chem. 24, 233 (1998). 

S. A. Mikhailov, Phys. Rev. B 65, 115312 (2002). 

M. Rontani, C. Cavazzoni, and G. Goldoni, Comp. Phys. 

Commun. 169, 430 (2005). 
™ M. B. Tavernier, E. Anisimovas, and F. M. Peeters, Phys. 

Rev. B 74, 125305 (2006). 
*° S. Tarucha, D. G. Austing, Y. Tokura, W. G. van der Wiel, 

and L. P. Kouwenhoven, Phys. Rev. Lett. 84, 2485 (2000). 



